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ABSTRACT 

The  classic  single  echelon  inventory  model  is  restruc- 
tured as  a  two  echelon  problem  in  which  demand  and  resupply 
are  deterministic.   Using  cost  minimization  as  the  objective, 
three  models  are  developed  which  address  the  problems  of 
(1)  no  stockouts  allowed  (EOQ) ,  (2)  backorders  allowed,  and 
(3)  finite  production  with  no  stockouts  allowed.   General 
solutions  for  the  optimal  policy  are  obtained  in  the  EOQ  and 
finite  production  models.   In  the  backorder  model,  the  ana- 
lytical argument  is  limited  to  the  case  in  which  only  time 
dependent  backorder  costs  occur.   Algorithms  are  developed 
for  solving  problems  for  all  three  models,  and  selected  para- 
meter values  are  used  to  test  the  behavior  of  the  models. 
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I.   INTRODUCTION 

Traditional  inventory  analysis  addresses  the  single 
echelon  inventory  problem.   In  this  problem,  customer  demands 
are  received  and  issues  are  made  from  a  single  outlet.   The 
outlet  in  turn  is  replenished  from  a  single  source. 

Unfortunately,  most  real  world  inventory  systems  are  not 
this  simple.   A  major  manufacturer,  for  example,  usually  has 
a  production  and  distribution  system  that  includes  (1)   fac- 
tories, (2)   factory  warehouses,  (3)   regional  warehouses, 
(4)   local  warehouses,  and  (5)   retail  outlets  [Hadley  and 
Whitin] .   Such  a  system  is  an  example  of  a  multiechelon  inven- 
tory system.   Each  level  of  the  system  is  a  separate,  distinct 
echelon.   Figure  1  depicts  graphically  this  organization.   The 
United  States  Navy  with  its  organization  of  inventory  control 
points,  supply  centers,  supply  depots,  shop  stores,  tenders  and 
underway  replenishment  ships  provides  an  example  of  a  large 
multiechelon  inventory  system. 

In  defending  their  single  echelon  models  Hadley  and  Whitin 
argue  that,  even  though  most  real  world  systems  are  multi- 
echelon, it  is  often  true  that  the  system  need  not  be  treated 
as  multiechelon.   They  contend  that  a  different  organization 
frequently  operates  each  level  (echelon)  of  the  system.   In 
the  example  of  Figure  1,  a  factory  and  its  associated  ware- 
houses might  be  operated  by  an  equipment  manufacturer  whereas 
the  regional  v;arehouses  might  be  separate  privately  owned 
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Figure  1.   Typical  Multiechelon  Inventory  System 


distributorships  not  under  the  administrative  control  of  the 
manufacturer.   In  like  manner,  the  local  outlets  might  be 
individually  owned  retail  stores.   Thus  we  would  have  a  sys- 
tem comprised  of  several  single  echelons  which  are  linked  by 
a  physical  dependence  but  separated  by  administrative 
controls. 

While  it  is  clear  that  there  are  many  examples  of  multi- 
echelon  systems  that  can  be  treated  as  a  series  of  independent 
single  echelon  problems,  it  should  be  equally  clear  that  there 
is  a  large  class  of  systems  which  must  be  treated  as  strongly 
dependent  multiechelon. 

It  is  the  objective  of  this  thesis  to  develop  and  investi- 
gate the  behavior  of  multiechelon  m.odels  constructed  from  the 
classic  single  echelon  inventory  models.   Specifically,  the 
following  deterministic  single-echelon  models  of  Hadley  and 
Whitin  will  be  restructured  as  two-echelon  models: 


a.  No  stockouts  allowed  (Economic  Order  Quantity) , 

b.  Backorders  allowed, 

c.  Finite  production  with  no  stockouts  allowed. 
Like  their  single  echelon  counterparts,  these  two-echelon 
models  will  be  single  item,  single  source  models. 


II.   THE  TWO  ECHELON  MODEL  WITH  NO  STOCKOUTS 

In  any  study  of  inventory  theory,  one  of  the  first  models 
investigated  is  the  deterministic  demand  model  with  no  stock- 
outs  allowed;  often  referred  to  as  the  Economic  Order  Quantity 
(EOQ)  model.   While  it  is  true  that  complete  deterministic 
demand  is  almost  never  known,  it  is  felt  that  the  mathematics 
of  this  model  provide  a  good  starting  point  for  any  inventory 
analysis.   Further,  it  is  felt  that  the  deterministic  model 
will  provide  an  insight  into  the  operation  of  a  stochastic 
demand  model.   It  is  considered  appropriate,  therefore,  to 
begin  this  multiechelon  analysis  by  considering  a  two  echelon 
extension  of  the  classic  EOQ  model. 

A.   MODEL  FORMULATION 

The  two  echelon  model  is  based  on  the  following 
assumptions : 

1.  The  upper  echelon  always  replenishes  its  supply  from 
the  same  outside  source.  The  lower  echelon  always  replenishes 
its  supply  from  the  upper  echelon. 

2.  The  upper  echelon  is  always  able  to  meet  the  demand 
of  the  lower  echelon. 

3.  External  customer  demand  always  occurs  at  the  lower 
echelon. 

4.  The  external  demand  rate  is  deterministic,  continuous, 
and  constant  with  a  value  of  X  units  per  year. 


5.  Neither  backorders  nor  lost  sales  are  allowed. 

6.  Procurement  lead  time  (PLT)  is  negligible, 

7.  The  order  quantity  at  the  upper  echelon,  Q2,  is  an 
integer  multiple  of  the  lower  echelon  order  quantity,  Ql. 
Let  n  =  Q^/Qt  /*  then  this  assumption  requires  n  >  i  and 
integer. 

The  optimal  policy  throughout  this  thesis  will  be  that 
which  minimizes  average  annual  variable  system  cost,  subject 
to  the  constraint (s)  of  the  model.   The  form  of  the  average 
annual  variable  system  cost  will  be  developed  by  first  deter- 

§ 

mining  the  system  cost  per  cycle,  and  then  dividing  this  cost 
by  the  cycle  length.   The  total  system  variable  cost  per  cycle 
is  the  sura  of  the  individual  echelon  variable  cycle  costs, 
K, (c)  and  K2 (c) .   The  cycle  length  is  defined  as  the  time  be- 
tween receipt  of  two  successive  orders  at  the  upper  echelon 
(echelon  two) . 

The  cost  per  cycle  at  echelon  two  is  comprised  of  an 
ordering  cost,  an  inventory  holding  cost  (IPIC)  and  a  purchase 
cost.  The  ordering  cost  is  assumed  to  be  a  constant  cost  per 
order  which  includes  the  administrative  costs  associated  v;ith 
inventory  review  and  order  (contract)  preparation.  This  cost 
is  independent  of  the  quantity  on  hand  or  on  order.  The  IHC 
includes  a  warehousing  cost,  an  obsolescence  cost  and  a  fore- 
gone opportunity  cost.  It  is  assumed  to  be  a  function  of  the 
inventory  on  hand  and  to  be  expressable  by 

T 


IHC   =   I  c/j  (t)dt. 
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Figure  2.   Two-Echelon  EOQ  Model. 

where  I-  is  the  average  cost  per  dollar  invested  in  inventory 
per  unit  time,  C  is  the  item  unit  cost,  T  is  the  cycle  length, 
and  c^^(t)  is  the  on  hand  inventory  at  echelon  two  at  time  t. 
The  purchase  cost  is  assumed  to  be  independent  of  the  quantity 
produced,  and  can  be  expressed  as  Q-C  where  Q^  is  the  quantity 
procured. 

The  variable  cost  per  cycle  at  echelon  two  can  now  be 
expressed  as 
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T 
K^ic)       =   A2  +  CQ2  +  I2C  J"c^2^t)dt  (2-1) 

The  height  of  the  on-hand  inventory  at  echelon  two  is 
Q2  -  Qi   =   [n-1]  Qi 

at  the  beginning  of  the  cycle  since  Q^  is  ordered,  arrives 
immediately,  and  is  partially  used  to  meet  an  immediate  demand 
of  Q,  units.   The  next  demand  at  echelon  two  is  for  an  amount 
Q,  which  occurs  Q,/X  time  units  later.   Similarly,  another 
demand  of  Q,  occurs  after  Q,/X  additional  time  units  have 
passed.   The  on-hand  level  at  echelon  two  continues  to  de- 
crease in  steps  of  Q,  until  it  reaches  zero.   It  remains  at 
zero  for  the  last  Q-,/x  time  units  of  T.   The  area  under  the 
on-hand  inventory  curve  for  one  cycle  of  echelon  two  is 
therefore 

-/o  m=l 

From  (2-1) ,  the  total  variable  cost  per  cycle  at  echelon 

two  is 

IjC   n[n-l]    Q^ 
K2(c)       =      A2    +   CQ2    +   2X •  ^^"^^ 

From  the  assumption  that  Q„  =  nQ,  ,  where  n  is  integer, 

it  follows  that  the  cycle  length  at  echelon  two  is  equal  to  n 

reorders  at  echelon  one.   Therefore,  from  Figure  2, 

Qo     nQ, 
^   =  A-  =   ^  •  ^2-3) 

The  average  annual  cost  at  echelon  two  is  obtained  by 
dividing  (2-2)  by  (2-3)  which  yields 
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Aa  I    C[n-1]Q^ 

K.       =      TT—  +    CX    + 


2  Q2  2Q2 

or,    as   a   function   of   Q,    and  n, 

AX                     I    CQ,  [n-1] 
^2      =      H^  -^   ^^    ^   2 •  ^2-4) 

Development  of  the  cost  per  cycle  at  echelon  one  follows 

an  argument  that  is  analogous  to  the  echelon  two  development, 

except  that  there  is  no  purchase  cost  incurred  at  echelon  one. 

Also,  from  the  assumption  that  Qp  is  an  integer  multiple  of 

Q,  ,  it  follows  that  in  each  cycle  there  are  n  reorders  at 

echelon  one,  i.e.,  K,(c)  is  linear  in  n.   Thus  the  cost  per 

cycle  can  be  written  as 

nI-,CQ^ 
Kj_(c)   =   nA^  +  — =^  .  (2-5) 

When  (2-5)  is  divided  by  (2-3)  the  average  annual  variable 
cost  is 

A,X    I  CQ^ 

The  total  average  annual  cost  of  the  system,  K,  is  the  sum 
of  (2-4)  and  (2-6)  . 

A,X    I,CQ,    Aa         I^C[n-l]Q, 

K   =   -i-  +  -i_-i  +  ^  +  CX  +  ~f . i  .  (2-7) 

Q^       2      nQj^  2 

The  model  has  assumed  that  the  unit  price  is  independent 
of  the  quantity  ordered.   Therefore  the  C  term  in  (2-7)  is  a 
constant,  and  can  be  dropped  from  (2-7)  resulting  in  (2-8) 
being  the  average  annual  variable  cost. 
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A,X    InCQ.    AX    i^C[n-l]Q, 
K   =   rri-  +  \         +  -4-  + 


Q- 


nQ. 


(2-8) 


B.   OPTIMAL  POLICY 

Determining  the  optimal  inventory  policy  involves  finding 
the  values  of  Q,  and  n  (call  them  Q*  and  n*)  which  minimize 
(2-8)  for  a  given  set  of  model  parameters.   To  find  this 
policy,  rewrite  (2-8)  as 

A, 


K(Q^,n)       =       [A^    +   ir]    ^Y  "^ 


[1^   +    I2[n-1]]CQ^ 


(2-9) 


Then  since  K(Q,,n)  is  convex  in  Q,  for  fixed  n,  Q?  (n)  must 


satisfy  the  equation 
9K 


A, 


3Q- 


or 


Qlin)      = 


2X  [A,    +   -=^] 
1        n    ■" 


[I,    +    I    [n-l]]C 
+   o =      0    r 


1/2 


[I^    +    I^ln-lJjC 


(2-10) 


After  substituting  (2-10)  for  Q,  in  (2-9)  and  collecting 
terms,  (2-9)  takes  the  following  form: 


K(Q*(n)  ,n)       =       {2X  [A^    +   ~]  [1^   +    I^ln-lJlC) 


1/2 


or 


K(Q*(n),n)       = 


2XA^I,C    + 


2XA2l-,C 
~n 


+      2XA,I^C[n-l]    +    2XA„I„C[^^] 

1     Z  Z     Z  Ti 


1/2 


Note  that  the  first  term  in  (2-11)  is  equivalent  to  K 


w 


(2-11) 


(Eq.  (2-11)  [Hadley  and  Whitin] ) .   In  fact,  it  is  readily  seen 
that  each  term  in  (2-11)  contains  a  form  of  this  classic 


formula.   Let 
.2 


K 


IJ 


2XA.  I  ,C  . 
1  3 
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Then  (2-11)  can  be  rewritten 

K  1/2 

K(n)  =  K(Q*  ,  n)  =  {K^  +  -—■  +   K^2t^-1^  +  ^22  ^^^  ^ 


(2-12) 


It  is  immediately  seen  that  (2-12)  is  positive  for  all 
feasible  values  of  n  (n  >.  1)  .   Further, 

lim  K(n)   =   [K^  +  K^,]""-^^  , 
n->l 

and  for  very  large  values  of  n,  ' 

1.    K(n)    ^   K^2  ^^'^^    >    0  ' 

8n        2n^/^ 

3     9^K(n)  ^   :^^12_    ^ 

:.  2        .  3/2       ^  • 
9n        4n 

If  the  assumption  that  n  is  an  integer  is  temporarily  re- 
laxed the  slope  of  K(n)  can  be  readily  investigated  by 

considering 

2      2 

21  ^   22  ^  ^2 

"2 ""2""     12 

8K   _nn  /oiox 

^  °     2K(n)1^2 •  (2-13) 

If  (2-13)  is  positive  when  n  =  1,  then  it  follows  that  it 
is  positive  over  all  n,  and  the  optimal  value  of  n  is  n*  =  1. 
If  (2-13)  is  negative  v/hen  n  =  1,  then  there  exist  values  of 
n  such  that  K(n)  is  less  than  K(l).   However,  since  it  has 
been  shown  that  (2-13)  is  positive  for  very  large  n,  it  follows 
that  n*  is  finite,  and  satisfies 
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8K 
9n 


=   0  , 


or 


n 


*  = 


"2      2  ' 

^21  ~  ^22 


K 


12 


1/2 


(2-14) 


A  further  consideration  of  the  case  where  (2-13)  is  nega- 
tive at  n  =  1  reveals  that 


2 


>   0 


over  the  range  of  n  values  from  1  to  beyond  n*.   Let  the  value 
of  n  where 


0 


be  denoted  by  n.   Then  in  the  range  1  £  n  <.  n  the  function 
K(n)  is  convex.   For  values  of  n  >  n,  (2-12)  is  not  convex. 
However,  (2-13)  will  always  be  positive.   Therefore  n*,  given 
by  (2-14),  is  the  optimal  non-integer  value  of  n.   The  reader 
will  note  that  the  above  arguments  indicate  that  (2-12)  is 
pseudo-convex . 

Since  K(n)  is  pseudo-convex  over  all  n,  it  follows  that 
K(n)  is  also  pseudo-convex  on  the  integer  values  of  n.   If  n* 
is  not  an  integer,  then  the  optimal  integer  value  of  n  will  be 
either  n  (the  smallest  value  greater  than  n*)  or  n  (the 
largest  value  less  than  n*) .   The  integer  optimal  will  be  the 
value  of  n  corresponding  to   K*  =  min  [K  (n )  ,  K  (n) ]  . 

In  summary,  the  optimal  inventory  policy  is  found  as 
follows : 
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2     2     2 

1.  Compute  K2-|,  ^22'  ^12  * 

2      2      2 

2.  Test  K2-,  -  K22  1.   K,2  •   If  this  inequality,  holds  then 

n  optimal  =1.   Go  to  step  8. 

2      2      2 

3.  If  Kp,  -  K22  *^  ^ip   then  compute  n*  using  (2-14). 

4.  If  n*  is  an  integer  then  n  optimal  =  n*.   Go  to 
step  8. 

5.  If  n*  is  not  an  integer,  compute  n  and  n. 

6.  Compute  K (n)  and  K(n)  from  (2-12). 

7.  K*  =  min  {K(n  ,  K(n)  }  and  n  optimal  is  the  value  of  n 
corresponding  to  K*  . 

8.  Compute  Q*  from  (2-10) . 

C.   NUMERICAL  EXAMPLE 

As  an  illustration  of  the  algorithm,  consider  the  follow- 
ing problem  where 

I^  =  0.75  A2  =  $200.00 

I2  =  0.50  X   =  100  units/year 

A^   =  $25.00  C   =  $100/unit 

2 
Begin  by  computing  K. .  =  2XA.I.C  , 


K'^^  =  3.0  X  10^ 
K22  =  2.0  X  10^ 
K^2  =  2.5  X  10^ 


Next  compare  K'^^   -  K22  =  1.0  x  10   >  K^^    =  2.5  x  10  , 
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so  compute  n*  using  (2-14).   The  value  is  2.0  which  is 
integer.   Therefore  n  optimal  =  2.   Using  (2-10)  the  optimal 
order  quantities  are 

Q*   =   166.67  units 

and 

Q*   =   166.67  X  2  =  333.34  units. 
2 

D.   PARAMETER  ANALYSIS 

The  primary  objective  of  this  analysis  was  to  observe  the 
effect  of  parameter  variations  on  optimal  n.   The  values  of 
four  of  the  five  parameters  were  fixed  and  the  value  of  the 
remaining  parameter  was  allowed  to  vary.   The  process  was  then 
repeated  for  each  of  the  other  parameters  so  that  the  effect 
of  each  could  be  observed. 

Representative  graphs  of  optimal  n  versus  a  given  para- 
meter are  illustrated  in  Figures  7  through  10.  In  Figure  7 
it  is  seen  that  n  optimal  is  inversely  related  to  A, .  How- 
ever, as  noted  in  Figure  8,  there  is  a  direct  relationship 
between  n  optimal  and  A„.  Further,  a  comparison  of  Figures 
7  and  8  reveals  that,  for  the  parameters  selected,  n  optimal 
is  much  more  sensitive  to  A„  than  to  A, . 

Figures  9  and  10  indicate  the  relationship  of  n  optimal 
to  I,  and  I„,  respectively.   The  results  indicate  that  n 
optimal  is  directly  related  to  I-.,  inversely  related  to  1.2' 
and  that  n  optimal  is  more  sensitive  to  a  change  in  holding 
cost  at  echelon  two  than  at  echelon  one. 


An  analysis  of  the  effect  of  demand  rate  on  n  optimal 
revealed  that  n  optimal  was  completely  insensitive  to  A  over 
the  range  5  <.  X  <.  10,000. 
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III.   FINITE  PRODUCTION  AT  THE  SECOND  ECHELON 

An  additional  two  echelon  variation  to  the  deterministic 
models  treated  [Hadley  and  Whitin]  addresses  the  problem  of 
finite  production.   This  model  would  have  application  in  any 
situation  where  the  upper  echelon  manufactures  as  well  as 
warehouses  a  given  product. 

A.   MODEL  FORMULATION 

The  model,  shown  graphically  in  Figure  3,  assumes  that 
echelon  two  produces  the  needed  material  at  a  constant  con- 
tinuous rate  of  ip  units  per  unit  time  {^    >.  X)  ,      The  cost  A„  is 
now  considered  to  be  a  set  up  cost  which  is  incurred  each  time 
a  new  production  run  is  initiated.   This  set  up  cost  is  inde- 
pendent of  the  quantity  produced.   The  rest  of  the  assumptions 
are  identical  to  those  of  the  EOQ  model. 

The  model  will,  as  in  Section  II,  seek  optimal  values  of 
Q-  and  Qy   which  minimize  the  average  annual  system  cost.   This 
cost  is  the  sum 

where  K.  is  the  average  annual  cost  at  the  ith  echelon.   In 
Section  II  it  was  shown  that 

^1    -    q7"  "^  i^i^'^i  • 

1 
To  develop  K2  it  is  necessary  to  define  the  cycle  length 


as 


T   =   tj_  +  T   +  T^  ,  V     .        (3-1) 
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where  t,  is  the  time  from  the  beginning  of  the  cycle  until  the 
start  of  production,  T   is  the  production  length,  and  T,  is 
the  time  required  to  deplete  the  inventory  at  echelon  two 
after  production  stops  (see  Figure  3) .   Note  that 

*i    -    X —  ~' 

P  "  X  ' 

and 


^d 


(Q2  -  Qi)    (Q2  -  Qi) 


De 


fine  i  as  the  number  of  times  that  the  quantity  Q,  is 


demanded  at  echelon  tv/o  during  the  production  period,  T  .   It 
can  be  seen  that 

t,  +  T  ^~' 
i   =   [    t   P]  (3-2) 

Ql 
where  t  =  t —  . 

A  \ 

As  in  Section  II,  the  variable  cost  per  cycle  at  echelon 
two,  K2(C),  will  have  the  form 

t,  +  T   +  T. 
1     p     d 

K2(C)   =   A2  +  I2C   A(t)  dt     .  (3-3) 

^1 

From  Figure  4  it  follows  that 
t,  +  T   +  T . 

r  1       p       ^ 

)^2^t)  ^t  =  f(^^  +  g^^  +  1)  +  h(i  +  2)  , 
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where  f(i)  is  the  area  under  the  first  i  saw  teeth,  g(i  +  1) 
is  the  shaded  area  shown  in  Figure  4  and  h(i  +  2)  is  the  area 
remaining  under  the  curve. 


Time 


Figure  4,   Finite  Production  Model,  Area  Under  the  Inventory 
Curve  at  Echelon  Two. 


From  Figure  4  it  can  be  shown  that 


Q^         [i-1]    qI^P        qI[^-X]  [i-1]  [i-2] 

f^i^     =     2^^  ZT--"  n 


(3-4) 


2X 


2\ 


and 


h(i)       = 


Q^    [n-i]  [n-i-1] 
2X 


(3-5) 
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From  Figure    5    the    area   g(i    +    1)    is    found    to   be 


g(i    +    1)       = 


-Q^    [n-1]^        qI^H^-I] 


2i/> 


2X 


Q^    [i[n-i]]         Q^{[i-1]  [ip-X]} 
+      _i +   —i 


(3-6) 


h"  ^  -H 


slope   ip/  ' 


K  b-H 
I 


Q- 


T 

d 


b      = 


~    [n-1]    -  j^    [i-1] 


c      = 


Ql    [i-l][^l 


d     =      ijjb 


Q^[n-1]    -   Q-^  ^    [i-1] 


Figure  5.   Finite  Production  Model.   Area  Under 
i  +  1^^  Saw  Tooth. 
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Substituting  (3-4),  (3-5)  and  (3-6)  into  (3-3)  yields 

I^CQ^    I„CQ^i|;[i-l] 


I2CQI  [n-i]  [n-i-1]         l2CQ^  [n-1] ^ 
■^  2X  2X 


I^CQ^ipli-l]  l2CQ^[i[n-i]] 


2X^  X 

I^C[i-l]  [4;-X] 

+  -^ 5 .  (3-7) 

X^ 

After  dividing  (3-7)  by  the  total  cycle  length  and  collecting 
terms  K2  can  be  written 

A2X         l2CQj^[n-l]  I2CQ^X  [n-2] 

K^   —   — ?: —  ■!■  o ~  !T~i •  (3  —  8) 

2      nQj^       2  24^ 

The  system  average  'annual  variable  cost  is  found  by  summing 
K,  and  K^  . 


X_ 
K(Qj_,n)       =       [^i+~]    +  — ^^ 2 ~ ~      ^^"^^ 


^2,  [Ii+l2[n-l]    -    l2^[n-2]]CQj^ 


B.   OPTIMAL  POLICY 

The  development  of  an  algorithm  for  determining  the  optimal 
inventory  policy  in  the  finite  production  model  follows  an 
argument  that  is  identical  to  the  one  discussed  in  Section  II. 
Since  (3-9)  is  convex  in  Q,  for  fixed  n,  Q|(n)  must  satisfy 
8K(Q, ,n) 
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or 


Q*(n)   = 


In 


1 


1/2 


^Cd^+I^tn-U  -  l2|[n-2) 


When  (3-10)  is  substituted  for  Q,  in  (3-9)  then 


K(n) 


^11  -^  ^+  ^12  ff^-1^  -  |t^-2n 


+ 


K^^[[n-ll  -  ^^[n-211 


n 


1/2 


(3-10) 


(3-11) 


where  KT  .  =  2XA.I  .C  . 
ID       1  3 

It  is  apparent  that  (3-11)  is  positive  for  all  values  of  n. 
Further, 

lim  K(n)   =   K?,  +  K?,  +  K?^^  +  ^l-X 
_^,  11     21     12\j;     221p  . 

and  for  very  large  values  of  n, 
1.   K(n)-^Kj^2""^'^^  >  0  ' 


2. 


3. 


dK 

^12 

8n    ' 

2nl/2 

a^K 

..    "^12 

9n     4n 


372 


>  0  , 


<  0  . 


If  the  assumption  that  n  is  an  integer  is  temporarily 
relaxed  the  slope  of  K(n)  can  be  investigated  using 


8K 

3n 


2  2 

21     2     X      22 
~T"  "^  ^12  f^"^^  "^  "T" 


2XK 


22 


n 


2  K(n) 


n 

ITT 


ipn 


(3-12) 
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If  (3-12)  is  positive  when  n  =  1,  it  is  also  positive  over 
all  n  and  it  follows  that  the  optimal  value  of  n  is  n*  =  1. 
If  (3-12)  is  negative  when  n  =  1,  then  there  exist  values  of 
n  such  that  K(n)  is  less  than  K(l).   However,  since  it  has 
been  shov/n  that  (3-12)  is  positive  for  very  large  n,  n*  must 
be  finite  and  satisfy 


or 


where 


and 


8K 
3n 


=   0  , 


n 


*   = 


2     ^2 

^21    ^22 


K 


12 


1/2 


(3-13) 


^2 

^22 


2  rij^-2X, 
^22^."!^ — ^  ' 


K 


12 


2  pip-X, 


A  further  consideration  of  the  case  where  (3-12)  is  nega- 
tive for  n  =  1  will  show  that 


8^K 


>    0 


over  the  range  of  n  values  from  1  to  beyond  n*.   Let  the  value 
of  n  where 


2 

a  K 

3n 


0 


be  denoted  by  n,  and  it  follows  that  in  the  range  1  <.  n  <.  n 
K(n)  is  convex.   For  values  of  n  >  n  K(n)  is  net  convex,  but 
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(3-12)  is  positive.   Therefore  n*  is  the  optimal  non  integer 
value  of  n  and  (3-11)  is  pseudo-convex. 

Since  (3-11)  is  pseudo-convex  over  all  n,  it  follows  that 
it  is  also  pseudo-convex  on  the  integer  values  of  n.   So  if 
n*  is  not  an  integer,  then  the  optimal  value  of  n  will  cor- 
respond to 

K*   =   min[K(n)  ,  K  (n)  ]  . 

The  optimal  policy  algorithm  can  now  be  written. 

2    ^^2        ^2 

1.  Compute  K^i /  ^2  2'  ^^^  ^1?  * 

2.  If  K^-j^  <.  K22  +  K^2  t^^^  ^*  =  1-   G°  to  step  8. 

2     ^2  '^2 

3.  If  K2-,  >  K22  "^  ^]  2  "'^^^^  compute  n*  using  (3-13). 

4.  If  n*  is  an  integer  then  n  optimal  =  n*.   Go  to  step 
8. 

5.  If  n*  is  not  an  integer  compute  n  and  n. 

6.  Compute  K(n)  and  K(n)  from  (3-11). 

7.  K*  =  min[K(n),  ^ (n) ]  and  n  optimal  is  the  value  of  n 
corresponding  to  K*. 

8.  Compute  Q*  from  (3-10). 

C.   SPECIAL  CASES  OF  i|; 
When  ^   =  \ 

A,X   IoCQt   A^X   I^CQ, 
K(n)  =  -J-+  2_1+  2  ^_J1 

Qq^     2      nQ^     2 
and 
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AX     [I  +1  ]CQ 

lim  K(n)   =  Q^  +   — ^^-j^ -'  (3-14) 

n-*-co  1 

Thus  it  is  seen  that  if  the  production  rate  is  equal  to  the 

demand  rate,  the  optimal  policy  is  to  start  the  production 

line  and  never  let  it  stop. 

As  one  would  expect,  when  ij;  -v  oo 

AX    I  CQ     AX 
lim  K(n)   =   ^r^  +  ^^    ^   +      ^ 


Qi       2      nQn 

l|;->-oo  ^1  ^1 


I  CQ, [n-1] 
+  — i 


and  it  is  immediately  seen  that  this  model  degenerates  to 
the  EOQ  model  of  Section  II. 

D.   NUMERICAL  EXAMPLE 

Consider  the  problem  where 

I    =   0.50  C   =   $100/unit 

I^   =   0.20  A   =   $50/unit 

X   =   25  units/year      A„  =   $200/unit 
t|/  =   75  units/year 

2    '^2  ^2 

First  compute  K.  .,  ^-y-y    ^"^^   ^i  ?  * 


21  ~   •"•"  ^  ■""  ^'11 

2  5  '^P 

22  =   2.0  K  10^  y\^ 

2  4  ^^2 

^2  =   5.0  X  10^  K^2 


K?^   =   5.0  X  10^  K?,   =   1.25  X  lo' 


K?^   =   2.0  X  10^  K?^   =   6.67  x  10^ 


2  4  ^^2  4 

k;^   =   5.0  X  10  K,^   =   3.33  x  10^ 
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Next  compare  K^^  with  K'^^  +  K^2~~^21  ""  ^'^  ^   ^^^  ^  ^22  "^  ^12  ^  "'■•^  ^  ^°^' 
so  compute  n*  using  (3-13).   The  value  is  3.64  which  is  not  an 
integer.   Therefore  n  =  4  and  n  =  3. 
From  (3-11) 

K(n)   =   $685.00, 

K(n)   =   $694.00, 
and  n  optimal  is  4.   Using  (3-10) 


and 


Q*   =   7.2  units 


Q*   =   28.8  units 


E.   PAFAMETER  ANALYSIS 

Figures  11  and  12  illustrate  the  effect  of  production  rate 
on  optimal  n.   If  the  holding  cost  at  echelon  two  is  small 
then  optimal  n  was  relatively  insensitive  to  a  change  in 
production  rate  (Figure  11) .   However,  the  results  of  Figure 
12  show  that  this  sensitivity  increases  as  1^    increases. 

The  responses  of  this  model  to  parameter  variations  were 
identical  to  the  responses  observed  in  Section  II. 
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IV.   THE  TWO  ECHELON  MODEL  WITH  BACKORDERS  PERMITTED 

A.   MODEL  FORMULATION 

The  EOQ  and  finite  production  models  were  based  on  the 
assumption  that  all  demands  would  be  immediately  satisfied. 
This  assumption  will  now  be  relaxed  in  order  to  investigate 
a  system  in  which  all  demands  are  ultimately  satisfied,  but 
where  it  is  permissible  to  accumulate  backorders  at  the  lower 
echelon.   No  backorders  are  allowed  in  echelon  two. 

In  this  model  dem^ands  v/hich  occur  when  the  lov/er  echelon 
is  out  of  stock  are  backordered  against  future  procurement. 
When  the  procurement  arrives  these  backordered  demands  are 
met,  and  any  excess  quantity  is  placed  in  stock.   Each  back- 
order  at  echelon  one  results  in  a  cost  of  the  form  tt  +  tt, 
where  tt  is  a  fixed  cha^rge  and  tt  is  a  time  dependent  cost. 

As  noted,  if  there  were  no  costs  associated  with  incurring 
backorders  [Hadley  and  Whitin] ,  then  it  would  be  optimal  to 
never  have  any  inventory  on  hand.   Conversely,  if  the  back- 
orders  are  sufficiently  expensive,  then  the  optimal  policy 
would  be  to  never  incur  any.   However,  for  an  intermediate 
range  of  backorder  costs,  it  will  usually  be  optimal  to  incur 
some  backorders  toward  the  end  of  the  cycle. 

With  the  exceptions  noted  in  the  preceding  paragraphs,  the 
model  is  predicated  on  the  same  assumptions  used  in  Section  II 

Define  s  as  the  number  of  backorders  at  echelon  one  when 
an  order  is  received.   Then  from  Fiqure  6  it  follows  that  the 
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on  hand  inventory  at  echelon  one  immediately  after  receipt  of 
an  order  is  (Q,  -  s)  units.   It  is  assumed  that  Q-,  >.  s.   Then 
the  echelon  one  on  hand  inventory  varies  from  (Q,  -  s)  to 
zero. 
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Figure  6.   Two-Echelon  Backorders  Allowed  Model. 


The  optimal  policy  continues  to  be  that  which  minimizes 
average  annual  variable  system  cost  subject  to  the  model 
constraints.   Since  the  model  does  not  permit  backorders  at 
the  second  echelon,  the  average  annual  variable  cost  at 
echelon  two  v;ill  be  identical  to  equation  (2-4)  ,    which  for 
convenience  is  restated  here. 
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A,X    I,CQ,  tn-1) 

Development  of  the  echelon  one  variable  costs  follow 
identically  the  single  echelon  backorders  permitted  model  of 
Hadley  and  Whitin.   In  the  interest  of  brevity,  this  develop- 
ment will  be  omitted  and  the  results  (equation  2-17)  [Hadley 
and  Whitin]  merely  stated. 

«1  =  57  ^  -^^ *  57  [-^^  \T-^  ■  '^-2' 

Summing  (4-I)  and  (4-2),  the  average  annual  variable  cost 


IS 


A^X  A^X     IiC[Q,-s]^  .  1   r  >   .  ^-'1 


l2CQj^[n-l] 


B.   OPTIMAL  POLICY 

The  optimal  solution  seeks  the  values  of  the  decision 
variables  Q*,  s*,  and  n*  which  minimize  (4-3)  subject  to  the 
constraint  that  all  demands  are  ultimately  satisfied.   To  do 
this,  note  that  (4-3)  can  be  rewritten  as 

^  =  I  ^  ^  -  ^x-  -  sf  -  lor  '        <^-^' 


where 


X   =   Ij^C  +  IT  , 

A 

A   =   A,  +  -^  , 
1    n   ' 


and 
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I   =   I^  +  I2  [n-l] 


Then  from 


^     =   0 


it  follows  that 


s      = 


I^CQ^ 


-    ttX 


X 


(4-5) 


or 


■nX    +   Xs 


I,C 


Next, 


9Ql 

-2AA     .     IC 

2         2 
2Q^^ 

27tXs 

results    in 

• 

°i  = 

— 

2AX    +    2ttXs 

+    Xs^ 

IC 

U. 

XsJ 

2q: 


=     0 


ni/2 


(4-6) 


(4-7) 


To  solve  for  Q*  and  s*  in  terms  of  n  and  the  system  parameters, 
substitute  (4-7)  into  (4-5)  and  (4-6).  After  collecting  terms, 
the  resulting  equations  are 

1/2 


Q*(n)   = 


2  2 

2A  X  -  77  X 

2   2 
ICX  -  I^  C 


(4-8) 


and 


s*(n)   = 


If. 
X 


2AXX  -  TT^X^ 


2   2 
ICX  -  1^      C 


1/2 


ttX 

x~  • 


(4-9) 


The  expressions  (4-8)  and  (4-9)  can  now  be  back  substituted 
for  Q,  and  s  in  (4-4)  to  yield  an  expression  for  the  average 
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annual  variable  cost,  K,  as  a  function  of  n  and  the  system 
parameters.   After  regrouping  the  terms  this  expression  can 
be  written 


[IC[I^C+7t]-IJC  ]  [2AX[I^C+Tr]-Tr^X    ]  I^CttX 

(4-10) 


K(n)  =   i i i +  ^^pr 

[I^C  +  7T]  -^l^^fT 


Equations  (4-8) ,  (4-9)  ,  and  (4-10)  are  valid  only  if 

2  2 
2AXIC  -  7T  X   >.  0 .   Otherwise  an  analysis  similar  to  that  of 

As 

Hadley  and  Whitin  is  required  for  the  case  when  tt  and  it  are 
positive.   Such  an  analysis  is  beyond  the  scope  of  this  thesis. 

Because  of  the  complex  form  of  K(n)  given  by  (4-10)  as  a 
function  of  n  no  attempt  was  made  to  evaluate  whether  or  not 
K(n)  was  convex  or  pseudo-convex  in  n  when  it  and  tt  are  both 
positive.   However,  an  example  is  presented  later  in  which  an 
optimal  value  of  n  >  1  is  obtained  for  it  and  tt  positive. 

If  IT  =  0  (4-10)  reduces  to 

I  c  1/2 

K(n)       =       [2AXIC   -    2AXIj^C   -^]         ,  (4-11) 

To   expand    (4-11) ,    let 

K?  .    =      2XA.I.C  (4-12) 

ID  1    D 


Then, 


2 

K      X 

K(n)       =   -472    K?1^   ^   K^^Xln-l]    +   -^i- 


K^^X[n-l]  ^  K^,    I,    C    -^/^ 

+   -=^ Kr^    It    C — 

n  11      1  n 
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( 


and   since   x  = 


I^C   +    TT,    K(n)    reduces    to 


K(n)       = 


'  2      '^  2      ^ 

-±^  +   -^ +   K^    [n-1] 

I    C+TT        n[I    C+tt] 


K^2l"-l' 


n 


1/2 


(4-13) 


Following  the  analysis  of  Section  II,  it  is  immediately 
seen  that  (4-13)  is  positive  for  all  feasible  values  of  n,  and 
that 


lim  K(n)   = 
n-^1 


2  ^     2    ^ 

I  C  +  7T 


1/2 


Also,  as  n  gets  very  large 

1.   K(n)->K^2  ^^'^    >  ^    ' 


2.   |2i^  J^  ^  Q 

8n   ^  1/2      ' 
2n  ' 


8^K(n)  ^  ^12   ^ 


8n' 


2n- 


Assume,  for  the  sake  of  argument,  that  n  is  continuous 
for  all  n  >  1.   Then 


•K 


21 


2 

8n2 


[I,C+7r]n 


2     ^22 

2  "^  ^12  ■•■  ~2~ 
n 


K(n) 


172- 


(4-14) 


and  it  is  obvious  that  for  n  =  1  (4-14)  is  non  negative  if  and 
only  if 
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I  c 
2      2      2        1 

^21  -  ^12  ^  ^22  f^  "^  "^^     ' 

77 


(4-15) 


If  (4-15)  holds,  then  (4-14)  is  positive  for  all  n.   This 

implies  that  over  all  n,  the  slope  of  K(n)  is  positive,  and 
n*  =  1. 

IF  (4-14)  is  negative  at  n  =  1,  then  there  exist  values 

of  n  such  that  K(n)  <  K(l).  However,  it  has  been  shov/n  that 
for  very  large  n,  (4-14)  is  positive,  and  it  follows  that  n* 
is  finite  and  satisfies 


or 


n 


B.    =  0 

9n      ^  ' 


*   = 


2    ^  2 

^21^  "  ^22 


[I^C  +  7T] 


^12  ^^1*^  +  ^J 


1/2 


(4-16) 


To  complete  the  argument,  note  that  when  (4-14)  is 
negative  at  n  =  1, 


2 

3  K 

3n2 


>  0 


over  the  range  of  n  values  from  1  to  beyond  n*. 
of  n  where 


Let  the  value 


2 

9  K 

9n 


=   0 


be  noted  by  n,  and  it  follov;s  that  K(n)  is  convex  over  the 
range  1  £  n  <.  n.   For  values  of  n  >  n,  (4-13)  is  not  convex, 
but  in  this  range  (4-14)  is  always  positive.   Therefore  n* 
given  by  (4-16)  is  the  optimal  non  integer  value  of  n.   The 
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reader  should  again  note  that  the  above  arguments  indicate 
that  (4-13)  is  pseudo-convex. 

Since  K(n)  is  pseudo-convex  over  all  n,  it  is  also  pseudo- 
convex  on  the  integer  values  of  n.   Therefore,  if  (4-16)  does 
not  yield  an  integer  value  then  the  optimal  value  of  n  can  be 
found  by  evaluating  K(n)  and  K(n)  as  in  Section  II. 

An  algorithm  for  finding  the  optimal  policy  (in  the  case 
where  tt  =  0)  can  now  be  stated. 

2     2         2 

1.  Compute  ^21'    ^7?'    ^^"^  ^12  * 

I  C 

2.  Test  K^^   £  [K^2  "^  ^22^  ^"^  "^  ~^^  *   ""-^  ^^^  inequality 

TT 

holds  then  n  optimal  is  1.   Go  to  step  8. 


2       2      2       "^1^ 

3.  If  K2-,  >  [K,2  +  ^22-'  ^ "^  "*"  "^^"^       then  compute  n*  using 

TT 

(4-16)  . 

4.  If  n*  is  an  integer  then  n  optimal  =  n*.   Go  to  step  ! 

5.  If  n*  is  not  an  integer  compute  n  and  n. 

6.  Compute  K (n)  and  K(n)  from  (4-13). 

7.  K*  =  min [K (n)  and  K (n) ]  and  n  optimal  is  the  value  of 
n  corresponding  to  K*. 

8.  Compute  Q?  from 


Qj(n)      = 


and    s*    from 


2X  [A,    +  -^] 
1        n 


C[Ij^l2C  [n-1]    +    Ij^TT    +    l2TT[n-l]] 


1/2 


s*(n)       = 


I   C  o* 

1^  ^1 


Ij^C     +     TT 


1/2 


H 

ZZ 

0.25 

H 

= 

$25.00 

*2 

= 

$100.00 

C.   NUMERICAL  EXAMPLE 

Consider  the  problem  where 

I^   =   0.75  7T   =   $5.00 

TT   =   $100.00/year 

X   =   100  units/year 

C  =   $100/unit 

2 

Begin  by  computing  K. .  =  2XA.I.C. 

K'^^      =   1.5  X  10^ 
K22   =   5.0  X  10^ 

K^2   =   1-25  X  10^ 

I  C 
2      2      2         1 
Next  compare  K^,  <  [K,^  +  ^o?-'  ^'^   "^  ""^r— ]  •   The  inequality  does 

77 

not  hold  so  compute  n*  using  (4-17).   The  value  is  1.69  which 
is  not  an  integer.   From  (4-13) 

K(n)   =   $735.00 
and 

K(n)   =   $620.00 

Therefore  optimal  n  =  1.   The  values  of  Qt    ,    Q^    ,    and  s* 
are  now  found  to  be 

Q*   =   12.9  units 

Q*   =   12.9  units 

s*   =   8.16  units. 
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D.   PARAMETER  ANALYSIS 

As  noted  previously,  computational  difficulties  preclude 
an  analytical  treatment  of  the  general  form  of  the.  backorder 
model.  However,  the  response  of  the  model  to  parameter 
variations  can  be  approximated.  This  can  be  accomplished  by 
considering  separately  the  two  functions  (3-1)  and  (3-2)  and 
ignoring,  for  now,  the  integer  restriction  on  n. 

Equation  (4-1)  can  be  rewritten 

A^X         l2^^^2  "  ^1^ 

and  since  the  second  term  is  never  negative,  it  is  obvious 
that  Kp  is  convex  in  Q^  .   Using  the  Hessian  it  can  be  shown 
that  (4-2)  is  convex  in  Q,  and  s  if  and  only  if 

2A^X  [I^C  +  ^]  >  [ttX]^  .  (4-18) 

Then 

A^X  I,C[Q^  -  Q,] 

+   r~-  +   -^ ^ i-  (4-19) 

^2 

is  convex  if  and  only  if  (4-18)  holds. 

It  can  be  shown  that  when  (4-19)  is  not  convex  the  optimal 

inventory  policy  will  be  realized  v;hen  backorders  are  not 

allowed  (s*  =  0).   Thus  if  (4-18)  does  not  hold,  the  problem 

can  be  solved  by  setting  s*  =  0  and  using  the  algorithm 

described  in  Section  II  to  get  Q?  and  n*. 
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For  the  case  where  (4-19)  is  convex,  optimal  Q, ,  Q^  and 
s  can  be  found  from 


8K 
BQ. 


^      =   0 
9s 


(4-20) 


Equation  (4-20)  will  be  satisfied  when 


Q|      = 


s*(n) 


2.  2 


2X  [nA,    +   A2]  [I-jC  +  tt]    -    niT  X 


n[I,C    +   7T    +    I^I^C    [n-1]    +    l2C7T[n-l] 


— \ 


1/2 


Ij^CQ*(n)  -  ttX 
I,C  +  77 


and 


(4-21) 


(4-22) 


Q*(n)   =   nQ*(n)  , 


(4-23) 


and  the  optimal  policy  can  be  found  from  the  following 
algorithm. 

1.  Set  n  =  1  . 

2.  Compute  Q?(n),  s*(n),  QS(n)  and  K(n)  using  equations 
(4-21)  through  (4-23)  and  (4-19). 

3.  Set  n  =  2  and  repeat  step  2. 

4.  If  K(l)  <  K(2),  stop;  n*  =  1  and  the  optimal  policy 
is  known. 

5.  If  K(2)  <  K(l),  set  n  =  3  and  repeat  step  2. 

6.  Continue  solving  K(n)  by  increasing  n  in  steps  of  1 
until  K(n+1)  >.  K(n).   Stop  as  soon  as  K(n)  <  K(n+1);  n*  = 
current  value  of  n  and  the  optimal  solution  has  been  found. 
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The  algorithm  described  in  the  preceding  paragraph  was 
used  to  observe  the  effect  of  parameter  variations  on  optimal 
n.   In  general,  the  results  were  similar  to  those  observed  in 
the  EOQ  model.   For  example,  n  optimal  varied  inversely  with 
A,  and  1^,    but  directly  with  K^   and  I, . 

The  sensitivity  of  n  to  the  various  parameters  was  highly 
dependent  on  the  selected  values  of  the  backorder  costs.   For 
example,  Figure  13  illustrates  the  relationship  between  optimal 
n  and  1^   in  a  situation  where  the  backorder  cost  is  small. 
Figure  14  illustrates  this  same  example  with  a  high  time  de- 
pendent backorder  cost.   The  difference  in  the  sensitivity  is 
obvious. 

Figures  15  and  16  illustrate  the  relationship  of  optimal  n 
to  a  change  in  A„ .   In  Figure  15  the  backorder  costs  were  low, 
and  optimal  n  was  insensitive  to  a  change  in  A^.   In  this 
example  optimal  n  remained  at  tv70,  and  s*  increased  as  A^ 
increased.   When  the  time  dependent  backorder  cost  was  high 
optimal  n  varied  directly  v;ith  A^  (Figure  16) .   However,  the 
change  in  optimal  n  was  not  nearly  as  large  as  in  the  EOQ  model 

Figures  17  through  20  illustrate  the  relative  magnitude  of 
Q*  ,  Q*  ,  and  s*  as  A^  and  I^  are  allowed  to  vary.  The  breaks 
in  these  curves  occur  at  points  where  optimal  n  changes. 
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V.   RECQMr4ENDATIQNS  FOR  FURTHER  STUDY 

The  analysis  of  the  backorders  permitted  model  should  be 
completed  for  all  values  of  the  backorder  costs  it  and  tt  in- 
eluding  the  case  where  tt  >  0  and  it  =   0.   Hadley  and  Whitin 
discuss  this  case  in  their  single  echelon  development  and 
conclude  that  if  tt  =  0   then  s*  is  either  0  or  infinite. 
Initial  investigations  indicated  that  in  the  two  echelon  model 
there  could  exist  a  finite  value  of  s*  >  0  when  tt  =   0.   How- 
ever, this  investigation  was  not  completed  and  no  conclusions 
were  reached. 

It  is  felt  that  the  analytical  argument  used  throughout 
the  paper  could  be  applied  to  the  general  backorder  case  where 
TT  and  TT  are  both  greater  than  zero.   However,  the  complexity  of 
the  equations  would  greatly  complicate  this  development. 

The  models  discussed  in  this  thesis  should  be  extended  to 
more  than  one  activity  at  each  echelon.   If  it  can  be  assumed 
that  all  activities  within  a  given  echelon  order  at  the  same 
time,  then  it  is  particularly  easy  to  include  multiple  activ- 
ities within  an  echelon.   For  example,  assume  that  there  are 
K  activities  in  echelon  one.   Then  in  the  EOQ  model  K,  would 


take  the  form 


k 

•^1  =    I 


i=l 


A, .  X  .    It  .CO, . 
li   1     li   li 


°ii 


If  the  values  of  the  parameters  A. ,  X  . ,  I, . ,  and  C  are  constant 

^  i'   i'   li' 

for  all  i,  i  =  l,...,k,  then  the  expression  for  K^  reduces  to 
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the  trivial  case  of 

K'   =   k  K^  , 

where  K.  is  found  in  Section  II.   The  expression  for  K^ 
remains  unchanged. 

The  problem  of  multiple  activities  at  the  top  echelon  is 
usually  not  a  very  interesting  one,  since  each  top  echelon 
activity  is  usually  responsible  for  supplying  specified  activ- 
ities at  the  next  lower  echelon.   Therefore,  unless  the  model 
provides  for  lateral  resupply  actions,  the  problem  of  multiple 
activities  at  the  top  echelon  can  be  reduced  to  the  sum  of  a 
series  of  independent  problems,  each  containing  one  top 
echelon  activity. 

A  related  problem  which  should  be  investigated  is  the  ex- 
tension of  the  models  to  more  than  2  echelons.   In  the  case 
where  each  echelon  is  limited  to  one  activity,  this  extension 
would  not  be  difficult.   Indeed,  given  the  assumption  that  there 
is  an  integer  ratio  between  the  quantity  ordered  by  successive 
echelons,  the  model  v/ould  probably  take  the  form 

m 
K    =      I    K 

:=i  J 

where  the  relationship  between  the  Q.'s  is 

D 

^2   =  ^1  •    Ql 

Q3   =   n^  •  Q2  =  n^  •  n^  .  Q^ 


J-1 

Q.  =  n .  ,  •  Q .  ,  =   n   n.Q,  , 
3  3-1    "^3-1    ^^^   1^1 

and  nK  is  an  integer  for  all  k,  k  =  l,...,j-l.   Because  of  the 
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integer  property  of  the  n,  ' s ,  the  solution  technique  outlined 

in  Chapter  II  could  be  used  to  solve  recursively  for  K. 

If  the  m-echelon  problem  allows  k.  activities  at  the  i 

^  1 

echelon,  i  =  l,...,m,  then  for  large  values  of  k.  and  m,  the 
sheer  magnitude  of  the  problem  would  make  its  solution  ex- 
tremely difficult.   In  fact,  it  is  doubtful  that  the  proce- 
dures suggested  in  this  paper  could  be  utilized  for  a  problem 
of  this  type.   However,  if  the  model  is- small  enough,  one 
could  at  least  get  a  "feel"  of  the  model's  behavior. 
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VI.   SUMMARY 

Three  of  the  deterministic  models  of  Hadley  and  Whitin 
have  been  restructured  as  two  echelon  models.   Equations  for 
the  average  annual  variable  cost  are  derived,  and  standard 
mathematical  programming  techniques  are  utilized  to  find  the 
optimal  inventory  policy.   The  optimal  policy  is  defined  as 
that  which  minimizes  average  annual  variable  cost  subject  to 
the  constraints  of  the  model. 

General  solutions  are  obtained  for  the  EOQ  and  finite 
production  models.   Because  of  the  complexity  of  the  cost  equa- 
tion the  analytical  solution  of  the  backorders  permitted  model 
is  limited  to  the  case  where  tt  =  0 .   However,  a  technique  is 
developed  which  can  be  used  to  find  the  optimal  policy  in  the 
general  case. 

A  parametric  analysis  is  conducted  for  each  model  in  order 
to  study  the  behavior  of  the  models  under  parameter  variations, 
The  behavior  of  the  curves  of  n  optimal  for  the  EOQ  and  finite 
production  models  were  identical  except  that  the  finite  produc- 
tion model  enjoyed  a  significant  cost  advantage.   The  behavior 
of  the  backorder  model  was  very  dependent  on  the  magnitude  of 
the  backorder  cost.   For  example,  if  the  model  had  a  low  time 
dependent  backorder  cost  then  parameter  variations  had  little 
effect  on  the  optimal  behavior  of  the  model.   However,  if  the 
backorder  cost  was  set  sufficiently  high  then  the  model's 
behavior  v/as  similar  to  that  of  the  EOQ  model.   Figures  7 
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through  20  illustrate  the  behavior  of  the  three  models  under 
varying  parameter  values. 
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Figure  17.   Q*  ,  Q*  and  s*  :  A2  Backorder  Model 
with  Lov;  Backorder  Costs. 
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Figure  18.   Q*  ,  Q*  and  s*  :  A„  Backorder  Model  with 
Pligh  Backorder  Costs. 
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Figure  19 


Q*  ,     Q*    and  s*  :  1^    Backorder  iModel  with  Low 
Backorder  Costs. 


60 


in 


^    \ 

>1 

V 

\ 

\ 

(0 

\ 

+J 

\ 

■t-\ 

\ 

c 

:i 

o 

O 

o 

• 

O 

o 

in 

rH 

rH 

II 

II 

II 

t= 

<    ^ 

4J 

w 

•H 

Q) 

C 

13 

(d 

in 

o 

O 

> 

r- 

in 

o 

in 

• 

• 

>-\ 

(N 

5^ 

o 

o 

<J> 

</> 

0) 

-P 

II 

II 

11 

II 

<1) 

g 

rH 

CN 

U 

r- 

td 

H 

M 

< 

^ 

ni 

fl4 

o 


61 


LIST  OF  REFERENCES 


1.   Hadley,  G. ,  and  Whitin,  T.  M. ,  Analysis  of  Inventory 
Systems ,    Prentice-Hall,  Inc.,  1963. 


62 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 

2.  Library,  Code  0212 
Naval  Postgraduate  School 
Monterey,  California  93940 

3.  Assoc.  Professor  A.  W.  McMasters,  Code  55  Mg 
Department  of  Operations  Research  and 

Administrative  Sciences 
Naval  Postgraduate  School 
Monterey,  California  93940 

4.  Assoc.  Professor  D.  A.  Schrady,  Code  55  So 
Department  of  Operations  Research  and 

Administrative  Sciences 
Naval  Postgraduate  School 
Monterey,  California  93940 

5.  LCDR  J.  V.  Dell,  S.C.  USN 
654  Landis  Street 

Chula  Vista,  California  92010 

6.  Chief  Naval  Personnel 
(PERS  IIB) 
Department  of  Navy 
Washington,  D.C.  20370 

7.  Department  of  Operations  Research  and 

Administrative  Sciences,  Code  55 
Naval  Postgraduate  School 
Monterey,  California  93940 


63 


Unclassified 


Secviritv  Cldssification 


DOCUMENT  CONTROL  DATA  -R&D 


(Security  classilicalion  ot  lllle,   body  of  abstract  and  indf^xing  ennotatior^  must  be  entered  when  Ihe  overall  report  (s  classified) 


1      ORIGINATING    ACTIVITY   (Corporate  author) 

Naval  Postgraduate  School 
Monterey,  California  93940 


2«.   REPORT    SECURITY    CLASSIFICATION 

Unclassified 


2b.    GROUP 


3   REPORT  TITLE 


A  DETERMINISTIC  MULTIECHELON  INVENTORY  MODEL 


4.    DESCRIPTIVE   NOTES  (Type  of  report  and.inctusive  dates) 

Master's  Thesis,  March  1972 


5     Au  TMOR'.SI  f Firs/ name,  rrxiddle  initial,  last  name) 


Jack  Vining  Dell 
Lieutenant  Commander,  Supply  Corps,  United  States  Navy 


6.    REPORT    DATE 


March    197  2 


7a.     TOTAL    NO.    OF    PAGES 


65 


76.    NO.    OF    REFS 


•a.    CONTRACT    OR    GRANT    NO. 


b.    PROJECT    NO. 


9a.    ORIGINATOR'S    REPORT    NUM8ER(S| 


9b.   OTHER  REPORT  NO(S)  (Any  other  numbera  ihat  may  ba  aaalgned 
Ihls  report) 


10.    DISTRIBUTION    STATEMENT 


Approved  for  public  release;  distribution  unlimited. 


11.    SUPPLEMENTARY    NOTES 


12.    SPONSORING   MILI  TARV    ACTIVITY 


Naval  Postgraduate  School 
Monterey,  California  9  39  40 


13.  ABSTRAC  T 


The  classic  single  echelon  inventory  model  is  restructured  as  a 
two  echelon  problem  in  which  demand  and  resupply  are  deterministic. 
Using  cost  minimization  as  the  objective,  three  models  are  developed 
which  address  the  problems  of  (1)  no  stockouts  allowed  (EOQ) , 
(2)  backorders  allowed,  and  (3)  finite  production  with  no  stockouts 
allowed.   General  solutions  for  the  optimal  policy  are  obtained  in  the 
EOQ  and  finite  production  models.   In  the  backorder  model,  the  ana- 
lytical argument  is  limited  to  the  case  in  which  only  time  dependent 
backorder  costs  occur.   Algorithms  are  developed  for  solving  problems 
for  all  three  models,  and  selected  parameter  values  are  used  to  test 
the  behavior  of  the  models. 


DD/„r..i473 

S/N    01 01 -807-681  1 


(PAGE    1) 


64 


Unclassified 


Security  Classification 


1-31408 


Unclassified 


Security  Classifiration 


KEY    wo  R  OS 


Inventory  Theory 
Two-Echelon 
Deterministic 
Optimization 


FORM 


,1473     <BACK) 


ROLE  WT 


LINK     C 


S/N     010) -b07-6S21 


65 


Unclassified 


Security  Classification 


A  -  3  1  ■!  0  9 


Th 
D2 
c. 


134015 


Thesis 

D2975     Del) 

c.l  A  determlnrstic  multJ 

echelon    inventory  model 


1     JUL7S 

-r     J  U  I.  7  5 

tomv  M 

7     FEB   90 


-      2336  1 


Thesis 
D2975  Dell 

c.l       A  deterministic  multi- 
echelon  inventory  model. 


134015 


'hesD2975 

A  deterministic  multiechelon 


inventory  m 


3  2768  001  02532  3 

DUDLEY  KNOX  LIBRARY 


